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ABSTRACT 



We use the SWIFT Burst Alert Telescope (BAT) sample of hard x-ray selected 
active galactic nuclei (AGN) with a median redshift of 0.03 and the 2MASS J 
and K band photometry to examine the correlation of hard x-ray emission to 

^ ! Eddington ratio as well as the relationship of the J and K band nuclear luminosity 

to the hard x-ray luminosity. The BAT sample is almost unbiased by the effects 
of obscuration and thus offers the first large unbiased sample for the examination 

Trt ', of correlations between different wavelength bands. We find that the near-IR 

nuclear J and K band luminosity is related to the BAT (14 - 195 keV) luminosity 
over a factor of 10^ in luminosity {Lj^ ?» L^^'^rp)aiid thus is unlikely to be due 
^ I to dust. We also find that the Eddington ratio is proportional to the x-ray 

Q\ [ luminosity. This new result should be a strong constraint on models of the 

^2 ' formation of the broad band continuum. 

O I Subject headings: galaxies: Seyfert, surveys, X-rays: galaxies 

OO 

O 



1. Introduction 

The estimate of black hole masses and the ratio of luminosity to the luminosity de- 
rived from the black hole mass (the Eddington ratio) are the fundamental parameters in 
the evolution of active galaxies. The Eddington ratio is defined as L/LEdd, where LEdd = 
(Mbh/Mq) X 1.3 X lO^^erg s~^ or the limit where gravitational and radiation pressure are bal- 
anced. It is important because it is a measurement of the ratio of energy output of a source 
to the physical limit of energy output. Understanding the energy output of active galactic 
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nuclei (AGN) is essential to understanding both the nature of AGN sources and the effect 
they have on their environment. Ther efore, measuring both Mnff and the Eddington ratio 
are the subject of much recent work (JKoUmeier et al\ l2006l : iNetzer fc TrakhtenbrotI 120071 : 
Woo &: Urryll2002l ). based primarily on optical samples and the com bination of x-ray an d IR 
sample s. These papers have either used data from the literature (IWoo fc Urrvl 120021). the 
SPSS (JNetzer fc TrakhtenbrotI 120071 ) . or deeper optical, x-ray and IR data (JKoUmeier et al. 



20061 ) to obtain conclusions about the relations hip of the lumino s ity to black hole mass and 
its evolution across cosmic time. For example, iKoUmeier et al\ ( l2006l ) conclude that there 
is strong evidence that supermassive BHs gain most of their mass while radiating close to 
the Eddington limit, and that the fueling rates in luminous AGNs are determined by BH 
self-regulation of the accretion flow rather than galactic-scale dynamical disturbances. 

Whereas previous samples are strongly biased against the most numerous AGN in the 



unive rse, those with high line of sight column densities, e.g. 'obscured' AGN (JMatt et al. 



20061 ). our BAT sample, since it is selected solely based on 14-195 keV flux, is not. However 
the BAT sample, as opposed to most optical samples, has a very low median redshift (0.03) 
and very few objects at z > 0.1. Si nce there is a strong relationship betwee n line of sight 
obscuration and intrinsic luminosity (JTozzi et a/.ll2006l : iMarkwardt et a/.ll2005l ) , such that the 
less luminous sources are much more likely to be obscured, the optical samples, which miss 
most of the obscured sources, are strongly deficient in low luminosity sources. Furthermore, 
there is often a complex relationship between measures of obscuration in the x-ray, IR, and 
optical bands and there may be stro ng evolutionary effects such that the fraction of obscured 
objects is higher at higher redshifts (JLa Franca et a/.ll2005l ). If there is a physical relationship 
between luminosity or stage in the life of the AG N and the occurrence of obscuration, as 
predicted by several models (JHopkins et aliUOOw . then the results from such surveys may 
be rather different than The BAT sample which is heavily weighted to the local universe 
with a mean redshift of ~ 0.03. 

In addition, there has been no conclusive determination of the origin of the 'non-thermal' 
continuum in AGN. In the x-ray band, there are two mechanisms that have received exten- 
sive discussion: thermal Comptonization of photons produced by the accr etion disk from 



a pop ulation of high energy particles, pe rhaps in a corona above the disk (IZdziarski et al 



1990l ). and emission due to jet processes (JKording. Falcke. fc Corbelll2006l ). However, there 
is at present no agreement as to the origin of 'non-thermal' IR and optical emission in radio 



quiet objects or even whether it exists. Most models 



have focused on disk models (e.g. ISun fc MalkanI (11989): ICzernv fc Elvis! (119871 )) with some 



Comptonization. Dusty torus models (jPier fc Krolik 



o r the origin of the opti cal emission 



19931 ), as the origin of the K band 



flux, have received support from reverberation observations (jSuganuma et al.l l2006l ) in the 
optical (UBV) and near-infrared (JHK) bands. These studies found a clear time-delayed 
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behavior of the K band with respect to the V band, strongly suggesting that the bulk of 
the K- band emiss ion originates from thermal emission of dust grains in an optically thick 
torus (IGlasd 120041) ■ There are strong correlations between the optical and x-ray propertie s 



( iBorosonl l2002l : iBrocksopp et al\ l2006l ) and the IR properties (JAlonso-Herrero et al 



200l|) 



whose origin is not easy to explain if these continuum components originate from completely 
different physic al mechanisms as is predicted by the most modern dusty torus models of the 
IR continuum (JNenkova et a/.ll2008l ). 



In this paper we present the first results on the Eddington ratio from an unbiased survey 
of AGN selected in the 14 - 195 keV band with the Swift BAT survey. We find a strong 
correlation of the Eddington ratio to the hard x-ray luminosity over a wide range in L/L^^^ 
(10~^ to 3) in our sample, as opposed to the narrow range in the sample of iKoUmeier et al 



( 120061 ) selected from the AGES survey. In addition we find that the luminosity of the point- 
like component in the 2MASS J and K band data is almost linearly correlated with the 
high energy luminosity indicating that these two components are very closely coupled and 
probably have the same physical origin. Since the x-ray emission is not due to dust processes, 
these results indicate most of the near IR is not either. 



2. Survey 



The BAT 9 month surve y (JTueller et a/.ll2008l ) is an all sky survey in the 14 - 195 keV 



band (JMarkwardt et a/.ll2005l ) with a sensitivity threshold that varies from 1— 3x 10 ^^ erg cm ^ s ^ 



in the 14 - 195 keV band, as a function of sky exposure. At \b\ > 15°, the survey is almost 
complete with 130 total sources of which 103 are identified AGN and the rest are galac- 
tic objects with only 1 unidentified source. Of the 103 AGN, 14 are Blazars and will not 
be included in the following discussion. All but 3 of the identified objects have redshifts. 
The survey is contin uing and the sensitivity is increasing roughly as the square root of time 
dTueller et allboosh . 



Black Hole Mass Estimation 



While there are several published papers on scaling bl ack hole masses to opti c al an d 
IR luminosities, we have chosen to use the recent paper by iNovak. Faber. fc Dekell (120061 ). 
They give scaling relations between the central black hole mass and global galactic properties 
based on optical and IR data. Since, for our sample, the only uniform optical or IR data 
set is the 2MASS data and we wish to minimize the contribution of possible AGN light to 
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the total IR light, we use J and K band magnitudes. As is well known, the ratio of the 
host galaxy to the QSO photon flux is at a maximum in the NIR, because the SED of many 
stellar populations has a maximum at 1.6 /^m and the AGN flux has a minimum at 1.2 /im 
( lElvis et a/.lll994j : [Richards et al\\200w . The AGN flux increases at shorter wavelengths due 
to accretion- disk or power-law emission and at longer wavelengths due to thermal emission 
of hot dust heated by the AGN and/or starformation contributions.. 

We obtained the total J and K 2MASS magnitudes (mj(tot) and mj^(tot)) for the BAT 
AGN sample, as listed in NED. Additionally, we used the IRAF task qphot to obtain circular 
aperture photometry of the central region, equal to the PSF FWHM of the 2MASS images. 
While we label these values as the nuclear magnitudes (mj(nuc) and mi^(nuc)), they are 
clearly overestimates to the nuclear flux, since no provision is made for the contribution 
of the galaxy light in the central region. We have found that in most objects the possible 
contribution to the total J band light from the nucleus is small and in most cases is well 
measured. This is important, since if the J band nuclear light dominated the emission, our 
estimates of BH mass would be severely impacted and the scatter in relations involving the 
BH mass (such as the Eddington ratio) would have much more noise. 

From the apparent magnitudes, we computed absolute magnitude s (Mj(tot), Mxitot), 
Mj{nuc), and Mxinuc)) with the redshifts listed in lTueller et al\ (120081 ) and assuming Hq = 
75kms~^ Mpc~^. From the total and nuclear absolute magnitudes, we computed estimates 
of the galactic/stellar component of the NIR emission {Mj{stellar) and Mxistellar)) as 
M,,.„„. = 2.51og(Y/(Y - 1)) + Mtot where x = Ft,t/F^^, = iQ-^AMtot-M^^.) _ Following 



l^.itellar 



Novak. Faber. fc Dekell (120061 ). we obtain black hole masses for each AGN in our sample 



using their relation: log{M{BH)) = 8.19 + 0.524 x (-Mx(.teHar) -23) or that logMf^g) 



1.3 log L{K stellar)- Our estimated black hole masses can be found in IWinter et all (120081 ). 
along with the 0.3-10 keV X-ray properties. 

W e note that aperture photometry in the J and K band (IWard et a/.lll988l : lGlass fc Moorwood 



19851 ) allows one to directly compare the AGN J,K band fluxes to the stellar flux for a sub- 
sample of our objects. We flnd that, in general, the total (stellar +nuclear) J and K band 
flux is dominated by star light with the AGN contribution, even for the most luminous 
BAT selected AGN, being less than ^ 30% of the total luminosity. This dominance of the 
near IR light by sta rs is beautifully conflrmed by the spectroscopic IR data for PG quasars 
(JDasyra et a/.l 120071 1. Even in these very luminous quasars (all of which have non-stellar 
IR luminosities more than 3 times larger than the stellar luminosities of the most luminous 
object in our sample and 80 times larg er than most of the objects), about 10% of the total H 
band light is due to stars. The work of JBentz et al\ (12006! ) . which contains some of the same 
objects as in our sample, shows that even in the optical band at 5100A (where the relative 
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dominance of AGN light is larger than in the near IR) none of their objects are dominated by 
AGN light. We thus believe that the noise introduced by the measurement error of the AGN 



light on the estimate of black hole masses for these objects is relatively small (jPeng et al 



20061 ). We have compared our photometry with that in the Peng et al. paper for both the 
extended emission and the nuclear emission. Peng et al. perform a detailed deconvolution 
of the 2MASS images and thus, presumably have done a more careful job than we have. 
However we find only ?« 1/3 mag scatter in the K band nuclear magnitudes between the 
Peng results and ours, indicating that our technique is sufficient to the task at hand and our 
relatively coarse photometry does not effect any of the results in this paper. 



Bolometric Correction 



We use a constant ratio of 15 to transfer the Swi ft BAT (14 - 195k eV) luminosity to 
bolometric luminosity, based on the broad band SEDS (JBarger et a/.ll2005l ) and the measure- 
ment of an E~^ photon index for the Swift BAT sources over the entire 14 - 195 keV band. 
While there seems to be a relat ion between the softer band x-ray luminosity and bolometric 
luminosity (ISteffen et adl2006l ). the effect is rather shallow and we choose to neglect it in 
this paper since it has not yet been measured in our harder x-ray band. If the effect is similar 
in the BAT band to that in the 0.5-2 keV band, it will tend to increase our estimate of the 
bolometric luminosity of the higher luminosity objects. We are accumulating broad band 
SEDS for all of these objects, but for many of the objects the nucleus is very weak in the 
UV and optical bands and thus measuring its luminosity is difficult. To obtai n good quality 



estimates of the nuclear optical and UV fluxes will probably require HST (JMalkan et al. 
19981 ) or adaptive optics data. 



5. Results 



5.1. Eddington Ratio 

Using the apparent magnitudes defined in Section 3, we derived associated NIR fluxes 
with the flux for zero-magnitude zero point conversion value s of 3.129 x 10~^° ergs~^ cm~^ 



for the J band and 4.283 x 10 ^^ erg s ^ cm ^ for the K^ band (jCohen et a/.ll2003l ). Luminosi 



ties, for example the nuclear luminosity Lj(„„c) and the galactic/stellar luminosity \^K{steUar)i 
were computed, as for the absolute magnitudes. Additionally, we computed the bolometric 
luminosity (bolometric correction xLbat) and the Eddington lumino s ity (u sing the esti- 
mated black hole mass from Mx (stellar) and the iNovak. Faber. &: Dekell (120061 ) equation) to 
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Fig. 1. — Plots of (a) the J filter nuclear luminosity versus the BAT x-ray luminosity (14 - 
195 keV), (b) the K filter stellar luminosity versus the BAT x-ray luminosity, (c) the ratio 
of the bolometric luminosi ty (estimated from the BAT l uminosity) to Eddington luminosity 
(estimated from M(BH) (JNovak. Faber. fc Dekell |2006| ) using the K filter stellar (total - 
nuclear) magnitude) (L/LEdd) versus the mass derived from the stellar K filter magnitude 
(in Mq), and (d) the ratio L/LEdd versus the BAT luminosity. The circles represent Syl - 
Syl.2 sources, the squares represent Syl.5-Syl.9, and the triangles represent Sy2s. For most 
of the plotted points, the error bars fit entirely within the size of the symbol. The lines in 
panels (a) and (d) represent linear fits to the data using the OLS bisector method (results 
are discussed in the text). 
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determine the Eddington ratio for each source. 

We show in Figure [T] (a) the correlation of BAT x-ray luminosity to Lj(n„c); (b) the 
relation of BAT x-ray luminosity to l-iK(steiiar), (c) the Eddington ratio versus black hole mass 
(M(BH)/M0), and (d) the correlation of the Eddington ratio to the BAT x-ray luminosity. Of 
the final 3 panels (b,c,d in Figured]), only the Eddington ratio to the x-ray luminosity shows a 
significant correlation, of the form logUo; /L^dd = (0.924±0.069) logLBAT+(-42.055±2.999) 
(solid line in Id) with a correlation coefficient of R2 = 0.45 for 98 points which has a 
probability of 3 x 10~^ by chance. Notice that both absorbed and unabsorbed AGN (primarily 
type II and I respectively) fall on the same scatter diagram indicating the absence of Compton 
thick objects, which mi ght be expected to have a high IR to x-ray ratio. Our results are 
similar to previous work (jWoo &: Urryll2002l ) where no correlation of Eddington ratio to black 
hole mas s is found. However, o ur median Eddington ratio of 0.069 is much less than that 
found by iKoUmeier et al\ (120061) who find valu es of ~ 0.2. We attribute this to the much 
higher luminosity of the iKollmeier et al\ (120061 ) sources, since we find for the most luminous 
objects in our sample a similar Eddington ratio of 0.2. The primary relationship is clearly 
with the BAT hard x-ray luminosity, since our assumption of a constant conversion from 
BAT to bolometric luminosity can only add noise, as does our method of estimating the 
black hole mass. 

The prime systematic errors are our estimate of the black hole mass and the bolometric 
correction. These should, on average, not reduce the significance of the relationship between 
the Eddington ratio and the hard x-ray luminosity. The sense of the errors is to have 
overestimated the black hole mass due to the improper correction for AGN light (e.g. the 
black hole mass is maximized in using the Nowak correlations if more of the K or J band 
light is assumed to be due to stars). This error will be a bigger problem for the higher 
luminosity objects because more of the K band light in the 2MASS aperture is due to the 
AGN. Since the hard x-ray luminosity does not have a similar possible error, we may have 
systematically underestimated the black hole mass and thus overestimated the Eddington 
ratio in the high luminosity objects. Similarly, if in the BAT band the x-ray light is a smaller 
fraction of the bolometr ic luminosity for lu minous objects, as it seems to be in the 0.5-2 keV 
band for type I objects (ISteffen et al\\200w . then the true bolometric luminosity is larger at 
the higher L(BAT). Both of these errors will tend to reinforce the sense of the correlation. 
Alternatively, if we have overestimated the AGN J band light it will reduce the total IR 
stellar light and thus our mass estimates and effect (see next section) the relationship of IR 
nuclear to x-ray light. Since this relation shows a scatter of only ^ 0.5 dex, this effect does 
not appear important. Of course there is some concern (as pointed out by the referee) that 
this correlation is an artifact of relating the x-ray luminosity to the x-ray luminosity divided 
by the black hole mass, where the black hole mass has a small range and the x-ray has a 



large range. We do not believe that this is the origin of our correlation since the correlation 
of L(BAT) with stellar mass has the same narrow range in mass and large range in L(x) 
and there is no correlation. Almost all correlations with black hole mass have the same bias, 
since there is only an approximately factor of 300 range in black hole masses estimated from 
bulge luminosities, one should exercise caution. 



5.2. Correlation of J and K band nuclei luminosities with hard x-ray 

As is obvious from Figure [H, the luminosity of the nucleus in the BAT and J band 
is strongly c orrelated, with a best fit (using t he ordinary least square bisector method as 
presented in llsobe. Feigelson. fc NelsonI (119861 )) of \og{Lj{nuc)) = logi^LsAT) x (1.257 ± 
0.077) + (—6.57 ± 3.09) and R^ = 0.63. The same correlation is seen with the K band 



nuclear luminosity (not shown). Using the same OLS bisector method. ISteffen et al\ (120061 ) 
finds a correlation between optical (2500A) and soft x-ray (2 keV) luminosity with a slope of 
0.721 ±0.011 in a sample of unobscured, optically-selected AGN. The non-linear relationship 
between optical and soft x-ray emission is partia lly true due to th e fact that low x-ray 
luminosity sources are more likely to be absorbed (ISteffen et a/.ll2003l ). Since the very hard 
x-ray/BAT and IR luminosities are relatively insensitive to obscuration, the almost linear 
relationship between these luminosities is not sensitive to obscuration effects and is thus 
more robust than the optical or UV to soft x-ray correlations. 



Discussion 



The median value of the Eddington ratio that we find is L/LEdd ~ 0.02 with a full 
range of ~ 1000. This valu e seern s to be at odds with the work of recent AGN samples, 
such as the iKoUmeier et all ( 120061 ) sample selected based on optical/IR/x-ray properties, 
which finds Lsoi/LEdd = 0.25 over a rather different luminosity range 10"^^ - lO^^ergs"^ 
in L(Bol) than the BAT sample. We believe that this is because the Kollmeier selection is 
really an optical one, since only the broad line sources have mass estimates. Su ch a sample is 
heavily biased towards "unobscured" sources which have higher luminosities (ISteffen et al. 
20031 ). Also, we now know that they also have higher L/LEdd (jWinter et a/.ll2008l ). Ad- 
ditionally, the Kollmeier selection is also biased towards much higher redshifts than the 
BAT survey, which could also explain the difference if, for instance, their is an evolution in 
Eddington rati os. The wide range of inferred Ed dington ratios from x-ray data has been 
noticed before (JMerloni. Heinz, fc di Matted 120051 ) and is also seen in other low z samples 
( Hopkins. Narayan. fc Hernquistll2005l ). However, our sample differs significantly from the 
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Ho. Filippenko. &: SargentI (119951 ) sample in not having objects with very low Eddington 
ratios. We are not sure why we would have a bias towards higher Eddington ratios than 
Ho. Filippenko. &: SargentI (119951 ). However, since we have not convolved our selection func- 
tion with a theoretical Eddington ratio model, one cannot be sure that our differences with 
Ho and Kollmeir are not due to selection biases, the luminosity range covered, and/or evo- 
lutionary effects. 



Softer x-ray surveys in the 2-10 keV band (e.g. ISteffen et al\ ( l2006l )) have found that 
the bolometric correction is a function of luminosity, such that the x-ray band carries less of 
the luminosity at higher luminosities. If this were also true for our sample it would increase 
the amplit ude of the correlation between hard x- ray and bolometric luminosities and change 
the slope (JHopkins. Narayan. &: HernquistI l2005l ) . 



The correlation between Eddingto n ratio and x-ray luminosity ha s , to ou r knowledge, not 
been reported from other samples (e.g. iMerloni. Heinz, fc di Matted (120051) ') but it is clearly 
present in them. If one graphs the lMerloni. Heinz, fc di Matted ( 2005) data one finds exactly 
the sa me correlation as in our plot. There are theoretical reasons (JHopkins. Narayan. fc Hernquist 
20051 ) to believe that the most luminous objects have the highest Eddington ratios and the 
lower luminosity objects have lower IjEdd- The Hopkins paper gives the Eddington ratio as 
a conditional probability distribution in which the width is a function of Eddington ratio - 
this is hard to interpret in our present framework and we have not tested the predictions of 
this model. 

We note that our distribution of ma s ses is weighted towards higher masses than pre- 
dicted by iHopkins. Narayan. fc HernquistI (120051 ) at low z, but we cannot conclude that this 
is a strong result since we have not convolved the iHopkins. Narayan. fc HernquistI (120051 ) 
prediction with our sensitivity limit. 

It is well known that there is a break in the luminosity function of AGN, the ori- 
gin of which is not clear. Below the break in the luminosity function, the median Ed- 
dington ratio from our sample is -2.1 and the variance is 0.35, while at higher L the 
median is -1.64 and the variance is 0.199. This is in agreement with the predictions of 



Hopkins. Narayan. &: HernquistI (120051 ). that the variance in the Eddington ratio should de- 
crease at higher L, and provides a possible explanation of the break in the luminosity function 
as a change in the mean Eddington ratio of the population. 

Our finding of a strong almo s t line ar correlation of the near IR luminosities with the 
BAT luminosity (see iHorst. et al\ (120071 ) for a similar result) argues that similar processes 
produce the near IR continuum and the hard x-rays. Since the mechanism producing the 
x-ray emission is unrelated to dust, these results do not favor dust reprocessing models for 
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the origin of the near IR continuum, unless the main driver of the dust reprocessing is simply 
and linearly connected to the x-ray continuum intensity. 

DHM acknowledges support from the National Aeronautics and Space Administration 
(NASA) under LTSA Grant NAG5-13102 issued through the Office of Space Science. 
This publication makes use of data products from the Two Micron All Sky Survey. 
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